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Abstract. We derive simple conditions for the stability or instability 
of the synchronized oscillation of a class of networks of coupled phase- 
oscillators, which includes many of the systems used in neural modelling. 



1. INTRODUCTION 

The general equations for a system of n coupled phase-oscillators are given 

by 

(1) =5,(01, 02,..., l<i<n, 

where 9i denotes the phase of the i-th oscillator, and the functions gi encap- 
sulate both the internal dynamics of the i-th oscillator and its coupling to the 
other oscillators. Since the 9i's denote phases, we assume that the functions 
Qi are 2tt periodic with respect to all variables. 

The study of the dynamics of networks of phase-oscillators is a formidable 
mathematical problem, and it is therefore important to isolate special struc- 
tures of such systems which facilitate the study of at least some aspects of 
their dynamics. Of course, in order for the work to be relevant to modelling, 
these structures must be natural in the sense that they reflect some mod- 
elling assumptions. In this work we identify a special structure for systems of 
the form ^ which makes the study of synchronization much more tractable. 
Moreover this structure is general enough to encompass many of the systems 
which have been used in neural modelling in the 'frequency domain' approach 

A synchronized oscillation of the system (^, is a solution for which 
e,{t) = e{t) for all 1 < i < n, with 

(2) 0{t + T) = e{t) + 27r Vt e R 

We shall say that the system is synchronized if there exists a stable synchro- 
nized oscillation (see section O for the definition of stability). The study of 
synchronization is an important part in the quest to understand the dynamics 
of networks of coupled oscillators ^2 ^| . 

We derive simple conditions for the stability/instability of synchronization 
for a class of systems. As we shall point out, these conditions generalize several 

1 



2 



GUY KATRIEL 



stability conditions derived in previous works for more restricted classes of 
systems. 

We first present a special case of the structure that we identify in this 
work, before turning to the general case. Assume that the functions gi have 
the special form: 

n 

(3) 5,(01, 0„) = h{e,) + c,,f{9„ Oj), 

Where h[a) is a 27r-periodic function describing the internal dynamics of the 
oscillators, /(a,/3) is a function which is 27r-periodic in both variables and 
describes the coupling of pairs of oscillators, and Cij are constants describing 
the strength of the influence of oscillator j on oscillator i. Let us note the 
modelling assumptions implied by such a form for the functions 

(i) The influence of all the oscillators on oscillator i is the sum of terms 
each one of which represents the influence of one of the other oscillators. 
This is known in the literature as the assumption of conventional synaptic 
connections. 

(ii) The type of coupling among pairs of oscillators, represented by the function 
/, is identical for each pair of oscillators, and only the strength Cij of the 
coupling is allowed to vary with i and j (including the possibility that Cij = 0, 
so that there is no influence of oscillator j on oscillator i). A more general 
model, still satisfying the assumption of conventional synaptic connections, 
would have / depending on the indices i, j, but such a model falls outside the 
scope of this study and does not admit the pleasing analytical results proved 
here. 

(iii) The internal dynamics of the oscillators, represented by the function h, 
are identical. 

Models of the form (O are already of great interest from a modelling per- 
spective but our results in fact apply to more general models 
in which the functions gi have the form: 

n 

(4) g^{el,...,e^) = s{9,,Yl ^j))- 

The function S{cj), y), 27r-periodic in the first variable, can be thought of as a 
cut-off function limiting the input to each oscillator (see |Z| |S] for such exam- 
ples), so its dependence on y (which wc do not restrict) might be sigmoidal or 
hump-shaped . We note that this structure relaxes the 'conventional synap- 
tic connections' assumption (i). Systems of the form are special cases of 
systems of the form (@J in which S{(j), y) = h{(j)) + y. 
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2. Existence of a synchronized oscillation 

We first deal with the rather trivial matter of existence of a synchronized 
oscillation for systems with the structure Q). 

Theorem 1. // there exists c such that 

n 

(5) Cij = c 1 < i < n 
and if 

(6) S{9,cf{9,9)) > V6I G K 

then there exists a synchronized oscillation 9i ~ 62 = ■■■ = On = 6 of 

n 

(7) e[ = s{e,,Y,<^vf{s^,^,)), i<t<n 

which is unique up to time-translations. 

Indeed, the synchronized oscillation is simply the solution of the scalar 
equation 

(8) 9' ^S(9,cf(9,9)). 

Condition guarantees that each of the equations reduce to ^ upon 
substituting 9i — 92 — ■■■ = 9n = 9, while condition 10) implies that each 
solution of satisfies |(2J for an appropriate T > 0. It is easy to see that (0 
is a necessary condition for the existence of a solution of JHll satisfying (0, 
hence for the existence of a synchronized oscillation. We note that the period 
T of the synchronized oscillation is given by the formula 

dO 

S{9,cf{9,9)) 

We make some remarks about the condition jSJ: 

(i) Condition says that the sum of strengths of the couplings of each 
oscillator to all the others does not depend on the oscillator, and it is easy to 
see that if this condition is violated a synchronized oscillation will usually not 
exist, because such a function would have to satisfy two independent scalar 
differential equations. 

(ii) Condition Q can be considered as an idealization - in a 'general' network 
we do not expect the sum of strengths of the connections to each oscillator 
to be exactly constant. As a consequence, in a 'general' network we do not 
expect a synchronized oscillation to exist at all. However, there are two very 
important classes of systems for which condition will naturally hold: sym- 
metric networks and regular networks, as we explain below. Moreover, we 
expect many natural 'large' networks to satisfy (|SJ) in an approximate way, 
leading to approximate synchronization - see below. 
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(iii) A symmetry of a network is a permutation r G Sn with the property 
that dj — Cr[i)T(j) for aU 1 < i, j < The set of symmetries of the network 
forms a subgroup G of Sn- It is easy to check that if the subgroup G is 
transitive, that is if for any 1 < i,j < n, there exists a t £ G such that 
t(?;) = j, then |SJ holds. Thus if we call a network symmetric whenever its 
symmetry group is transitive, we have that any symmetric network satisfies 

A simple case is a cyclically-symmetric ring of oscillators, in which the 
symmetry group is generated by the permutation r(i) ^ i + 1 {mod n). For 
more on symmetry in networks of oscillators see jB]- 

(iv) We note the trivial fact that a linear combination of two matrices satisfy- 
ing ISJ also satisfies (0), so that our condition holds, for example, if we have 
two independent sets of connections, each of them with a different transitive 
symmetry group. 

(v) Consider the case of a network in which the connections have either 
strength s or strength (i.e., no connection). Such a network can be rep- 
resented by a directed graph in which the oscillators are the nodes and an 
arrow from node j to node i indicates that Cy = s. If this graph is k-regular 
in the sense that each node has k incoming arrows, then condition Q is sat- 
isfied. This includes 'random' networks generated by asigning k oscillators to 
each oscillator. 

(vi) If we relax the assumption of regularity and generate a random network 
by a (perhaps more natural) process of setting Cy = s with probability p and 
Cij = with probability 1 — p, we can still expect that when the number n of 
oscillators is large, the network will satisfy © in an approximate way. Indeed 
by the central limit theorem the sums in ^ will be approximately pns, with 
an error of 0(--^). Thus such a network can be considered as a small pertur- 
bation of a network satisfying ij^l with c = pns. Under such circumstances, we 
expect the perturbation to change the synchronized oscillation into a nearby 
entrained oscillation, that is a solution satifying 9i{t + T) — Oi{t) + 2tt for all 
1 < i < n. Moreover we expect the functions 9i to be close to one another, 
though not identical. We thus have the phenomenon of 'near-synchronization', 
and we also expect that the stability of the 'near-synchronized' oscillation will 
be the same as that of the synchronized solution of the system satisfying |3) 
off which we perturb. These heuristic remarks need to be backed by a careful 
analysis, but we leave this to future work. 

We formulate the corollary of theorem^ applied to systems of the form (01: 
Theorem 2. Assume 0) holds and 



(10) 



h{e) + cf{0,0) > e M 
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then there exists a synchronized oscillation 9i ^ 62 = ■■■ — 9n — Q of 

n 

(11) 0,; = l<i<n 

which is unique up to time-translations. 



3. STABILITY OF THE SYNCHRONIZED OSCILLATION 

We now turn to the more interesting issue of stability of the synchronized 
oscillation. 

We recall some standard definitions. Let d denote the euclidean distance in 
R". A solution 9{t) = {Oi{t), ...,9n{t)) of (QJ is said to be Liapunov stable 
if for any e > there is a i5 > such that whenever to G I^- ^'o G I^" and 
d{9o, 9{to)) < S, we have d{9{t), 9{t)) < e for aU t > 0, where 9{t) denotes the 
solution of with initial condition 9{0) = 9q. 

9{t) is said to be (asymptotically) stable if it is Liapunov stable and there 
exists a (5i > such that whenever tp G M, G K" and d{9o,9{to)) < Si, 
there exists some <i £ K such that limt^oc d{9{t) , 9{t — ti)) = 0, where 9{t) 
denotes the solution of ^ with initial condition 9{0) — 9q. 

9{t) is said to be unstable if it is not Liapunov stable. 

We note that we are dealing here with local stability, so that the fact that 
the synchronized oscillation is stable does not imply that synchronization 
will be achieved from any initial state, but only that there is some open set 
of initial conditions that will lead to synchronization. The study of global 
synchronization properties of systems of phase-oscillators is a very intriguing 
problem. 

In general the stability of a synchronized oscillation is to be determined by 
examining the time-periodic linear system obtained by linearizing ^ around 
the synchronized solution: 

(12) u' = A{t)u, 1 < i < n, 
where Ait) is the n x n matrix with 

(13) A,it) = ^^m,-,m), 

which in general can only be done by numerical computation. The crucial 
point is that the special structure I^J entails a special structure for the lin- 
earized system (|T^ , which makes it possible to proceed much further with the 
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analytic investigation than in the general case. Differentiating we have 
% 



(14) 1^(01,.. .,0„) = 



Thus substituting 9i = 02 = ■■■ = On = 6 ((1^ wc obtain 

(15) A{t)^a(e{t))I + bie{t))C, 

where / is the identity matrix, C is the connection matrix Cij = cij , and the 
scalar functions a{9),b{9) are given by: 

(16) ai9) = ^{9, of {9, 9)) + c^{9, of {9, 9))^{9, 9) 

(17) b{9)^—{9,cf{9,9))^{9,9). 

Let us assume for the moment that the connection matrix C is diagonaliz- 
able: 

(18) C = M-^DM 
with D the diagonal matrix 

D = diag[Xi, A„], 

where , 1 < k < n are the eigenvalues of C. We note that Q implies that 
(1, 1, 1) is an eigenvector of C with eigenvalue c, so that we may assume 

(19) Ai = c. 

Let F{t) denote the nxn matrix- valued function which is the fundamental 
solution of (22), so that F{0) ^ I and 

(20) F' = A{t)F. 

We recall from the stability theory for periodic solutions of autonomous sys- 
tems that the synchronized oscillation is stable if all the eigenvalues of F(T) 
(T denoting the period of 9), known as the characteristic multipliers, are in the 
interior of the unit disk in the complex plane, except for one eigenvalue which 
is 1 (related to time-translation invariance). If at least one of the eigenvalues 
is in the exterior of the unit disk, the synchronized oscillation is unstable (see, 
e.g., [12], Ch. V, theorem 8.4). 
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Defining G = MFM-^, wc have from (HSJ, ^ 

(21) G' = [a(e{t))I + b{e{t))D]G, 

with G(0) ~ I. Since the linear system (|21(l is diagonal, its solution is 

Git) ^diag[n{t),...,r„it)] 
where the functions satisfy the scalar equations 

(22) r',{t)^[am) + Kbm)]nit), 
with ri(0) = 1, whose solution is 

r,{t) = e^o [a(eis))+\Meis))]ds ^ 

In particular we have 

G{T) = diag[e^o la(e{t))+Xib(e{t))]dt ^ ^^^^ [a(e(t))+A„b(e(t))]dt]^ 

so that the eigenvalues of i^(r) are 

(23) [Li^e''\ l<i<n. 
where 

(24) a, = / [a(eit)) + \,b(eit))]dt 



By making the change of variable 9 = 9{t), so that, using ©, d9 
S{9,cf{9,9))dt, we obtain 

In particular, for i = 1, using (|19|l . 

and using the easily-verified fact that 



±[S{9,cf{9,9)]^a{9)+cb{9), 



we obtain 



d d9 f^"" d 



-^=X ^[^(^--/(^-^))] 5(.,c/(M)) -X 
Thus /ii = 1 is the eigenvalue of F{T) corresponding to time-translations, and 
the sufficient condition for stability of the synchronized solution is that all the 
other eigenvalues, given by (|23|) . are inside the unit disc. This is equivalent 
to the condition that the numbers ai have negative real parts for 2 < i < n. 
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We can express ai in a somewhat simpler form by noting that from H26|l 
and the fact that cti = we have 



2tt 



a{0)d0 f^"" b{0)d0 



lo Si9,cfi9,9)) Jo Si9,cfi9,9)) 
so we can rewrite (|25|l as 

_ f^^ b{9)d9 

^'-^^^ 'Vo S{0,cf{0,0)y 

RecaUing the definition (|17() of b{9), we finaUy obtain the expression: 

(27) a.,^(A.-c) r ^i9,cfi0,0)&0,9) f 1 < . < n. 

Jo oy d(3 S{0,cf{0,0)) 

We note now that ahhough H27|) was derived under the assumption that 
C is diagonahzable, it continues to be vahd for any matrix C satisfying 
because the set of diagonahzable matrices is dense and both sides of (|27|l are 
continuous in C. We thus have: 

Theorem 3. Assume |3|) and @) are satisfied. Let 

(28) x = x{SJ,c)= ^(9,cf{9,9))^{9,9)- 



dy' ' ' "dp' ' ' Si0,cfi0,0))- 

(29) x[Re{Xi) - c] < 0, 2<i<n 

then the synchronized oscillation of is stable. 

(a) If for some 2 < i < n, x[i?e(Ai) — c] > then the synchronized oscillation 
of is unstable. 

We formulate the stability result for the special case of conventional synap- 
tic coupling: 

Theorem 4. Assume and i) i^) are satisfied. Let 



^^^^ ^ -v-,.,--y Qp^-^"'h{0) + cf{0,0y 

(i) If 1^29\) holds then the synchronized oscillation of 1^11]) is stable. 

(a) If for some 2 < i < n, x[i?e(Ai) — c] > then the synchronized oscillation 

of J^ll\) is unstable. 

The calculations involved in determining the stability of the synchronized 
oscillation using theorems|nior01split into two parts: calculation of the integral 
X, and study of the spectrum of the connection matrix C . The following 
reformulation of theorem |31 highlights the independence of these two aspects: 
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Theorem 5. Assume hold, and let x be defined by J^) . 

(31) min Re{\j) < c < max ReiXA 

2<i<n 2<i<n 

then the synchronized solution of any system of the form ^ with x 7^ is 

unstable. 

(n) If 

(32) min i?e(A,) > c 

2<i<n 

then the synchronized solution of is stable if X ^ ^ unstable if ^- 

(ill) If 

(33) max Re{Xi) < c 

2<i<n 

then the synchronized solution of ^ is stable if x > o,nd unstable if x < 0. 

We note that x depends on the connection matrix only through the 'overall 
strength' c. 

Let us study, as an example, the system 

n 

(34) 0',=uj- sm{0,) cos(0,), 

considered in ^U]. This is a system of the type with 

h{a) = uj, f{a, (i) ~ — sin(a) cos(/3). 

We assume Q holds. In order to satisfy H1U|I we must have w— c sin(0) cos(6') > 
for all 9, which is easily seen to be equivalent to the requirement 

(35) > M. 

Thus if (jSJ and I35|l hold, there exists a synchronized oscillation of 134|l com- 
puting X according to H3U|I we obtain 

■ 2 //IN 

sm (a) — -— -, 

Q Lo — csin(6') cos(6') 

and we see that x > whenever (|35|l holds. Thus from theorem 0] we obtain 
that the synchronized solution is stable iff the real parts of the eigenvalues of 
C (other than c itself) are all smaller than c, or, in other words, when all the 
eigenvalues oi C — cl have negative real parts. 
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4. Connection matrices and the stability of synchronized 

oscillations 

The spectrum of the connection matrix C is, of course, totaUy independent 
of the nonhnearities ft., /, S. This fact allows us, for some connection matrices, 
to conclude that the synchronized solution is unstable, independently of the 
internal dynamics of the oscillators and of the coupling functions, as shown 
by part (i) of theorem 

In this section we consider the role of the connection matrix through some 
examples, using some well-known results from linear algebra. 

4.1. Two-oscillator networks. In the case rt = 2 it is easy to see that if 
we assume no self-coupling, (cn = C22 = 0) the general system Q is in fact 
no more general than the particular case Hll|l . and in this case we may also 
absorb h into the /, so we are dealing with the system: 

(36) e[ = f {6^,92) 

The connection matrix is 



C 



1 

1 



and we have Ai = c = 1,A2 = —1 Thus H33() holds and from theorem [S] we 
obtain: 

Theorem 6. For a system of two coupled phase-oscillators of the form ^clb'}) 
for which 

f{e,e) > V6' e M, 

let: 

'■2'^ df ,^ ^. dd 



^ ' dp^-'"'f{9,9)- 
The synchronized solution is stable if X > ^ '^^'^ unstable ifx<0. 

4.2. Three-oscillator networks. We assume that there is no self-coupling, 
so that the connection matrix can be written in the form 

(37) C = 

By computing the characteristic polynomial of C'—cI, we find that the product 
of its two eigenvalues (besides the eigenvalue 0) is 

(A2-c)(A3-c) = A 

where 

(38) A = 3c^ - (ci2 + C23 + C3i)c + (C12C23 + C12C31 + C23C31), 
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and that their sum is 

(A2 - c) + (A3 - c) = -3c. 

Assume A < 0. Then A2,A3 must be real (because otherwise they are 
conjugate which would imply that the product (A2 — c)(A3— c) is non- negative), 
and (A2 — c), (A3 — c) are of opposite signs, hence part (i) of theoremElimplies 
that the synchronized oscillation is unstable. 

Now assume A > 0. If A2,A3 are non-real, hence conjugate, then we have 

3 

Re{X2) — c = Re{\z) — c ~ ^2'^ 

so that H32|l holds if c < and 1)33(1 holds if c > 0. If A2, A3 are real then 
A > implies that (A2 — c), (A3 — c) are of the same sign, which is opposite 
to the sign of c, so that again ((32|l holds if c < and H33|l holds if c > 0. 
We thus obtain: 

Theorem 7. Consider the system ^ with n ~ 3 and connection matrix given 
by jST)) , and assume that 0) holds. Define x by i2f^) . 

(i) If A. < and X (where A is given by 1^3 yp } then the synchronized 
oscillation is unstable. 

(a) If A > then the synchronized oscillation is stable if cx > and unstable 
if cx < 0. 

4.3. Larger networks. In principle, for larger networks, with theorem [Sj 
one can obtain algebraic conditions involving the elements of the connection 
matrix and x for determining the stability/instability of the synchronized so- 
lution, by applying the Routh-Hurwitz conditions to the characteristic poly- 
nomial of C — cl. Obviously these conditions become more complicated as the 
size of the network gets larger, and we shall not proceed along these lines but 
rather consider some specific classes of connection matrices for which simple 
results are possible. 

4.4. Non-negative connection matrices. A very important case is that in 
which all the entries of the connection matrix are non-negative. In this case 
we can apply the Perron- Frobenius theorem to obtain stability results. In 
fact to apply the Perron-Frobenius theorem we need to assume also that the 
matrix C is irreducible, which means that there does not exist a partition of 
the indices {1, n} into two disjoint sets A, B such that i € A, j G B implies 
Cij = 0. In other words, we assume that we cannot partition the oscillators 
into disjoint sets A, B so that oscillators in B have no influence on oscillators 
in A (we note that a simple sufficient condition for C to be irreducible is that 
Cij ^ for all i j). The Perron-Frobenius theorem states that if C is non- 
negative and irreducible then it has a unique eigenvalue A which is positive 
and has an eigenvector whose components are positive. Moreover A is a simple 
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eigenvalue and all other eigenvalues of C are less than or equal in modulus to 



Since in our case we already know of the eigenvalue Ai = c which has an 
eigenvector (1,1,..., 1) with positive components, the uniqueness part of the 
Perron- Frobenius theorem implies that all other eigenvalues satisfy Xi ^ c and 
I Ail < c. These two facts imply that Re{Xi) < c, so that holds. We thus 
conclude: 

Theorem 8. // the connection matrix C satisfies |^ and is non-negative and 
irreducible, and ^ holds, then the synchronized solution of is stable if 
X > and unstable ifx<0, where x is defined by 

Of course an analogous theorem is true for the case of non-positive matrices, 
where the roles of the conditions x > 0, x < ^re reversed. 

A special case of theorem |H1 was proved in ;4 (section 2.5), with x given 
in a somewhat less explicit form. There it was assumed that the system is 
of the form lll|l with h a constant and that / is of the special product form 
f{a,(3) = fi{a)f2{P)- The unneccesary assumption that C is a symmetric 
matrix was made. 

We note that theorem [S] includes the case of global uniform coupling given 
by Cij = s for all i,j (or for all i ^ j). This case has been treated, for 
the special case of (|ll|l in which the coupling f{a, (3) depends only /3, in 
0, and for the special case in which the coupling is of the form f{a,/3) = 
fo{(3 — a) + fi{a)f2{l3) and h is constant, in 2 . The stability results for 
synchronized oscillations in these works are corollaries of theorem [S] 

4.5. Cyclically symmetric systems. We now consider the case of a ring 
of oscillators with cyclic symmetry, which means that the effect of the j-th 
oscillator on the j-th oscillator is the same as that of the i + k-th oscillator 
j-|-fc-th oscillator (where i+k, j+k are considered modulo n). In mathematical 
terms this means that the connection matrix C is circulant that is, one 
having the form 



with d G some vector. As we noted before, a circulant matrix automati- 
cally satisfies (jSJ. 



We have the following formula for the eigenvalues of a circulant matrix, 
where p — e^^'^*/": 



A. 



(39) 



'ij ^(j — 2+1) 7nod 



n 



(40) 




m— 1 



This implies 



(41) 




m— 1 
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The formula (|41|l . together with the identity 1 — cos(x) = 2sin^(|), give 
us the foUowing stabihty resuh for cychcaUy symmetric matrices: 

Theorem 9. Assume ^ is satisfied, and assume that the network is cyclically 
symmetric, its connection matrix given by Let x be defined by and 

define, for 1 < k < n: 

d„ sm^ [tt^ '-^ '-j . 

m—l 

(i) If some of the 's are postitive and some are negative, the synchronized 
oscillation of is unstable. 

(a) If Xlk > for all 1 < k < n then the synchronized oscillation is stable. 
(Hi) If Xlk < for all I < k < n then the synchronized oscillation is unstable. 

4.6. Two populations of interacting oscillators. Let us consider 2n oscil- 
lators, subdivided into two groups, each of size n, with the pairs of oscillators 
within each group coupled with strength p, and any pair of oscillators belong- 
ing to different groups are coupled with strength q. 



Cij — p if 1 < i ^ j < n or n + 1 < i ^ j < 2n 
= q if ^ l£ "i ^ n + I < j < 2n, 
= q if n + 1 < i < 2n, 1 < j < n, 
= if i = j- 
For example in the case n — 3 the connection matrix is 

/O p p q q q\ 

p p q q q 
^^ppQqqq 

q q q p p 

q q q p p 
\q q q p p 0/ 

Such a matrix satisfies (O with c = {n—l)p+nq. Let us denote by G M^" 
(1 < A; < 2n) vector with 1 in the fc-th place and elesewhere. It is is easy to 
verify that the eigenvalues of C are c with eigenvector (1, 1, 1..., 1), {n—l)p—nq 
with eigenvector 1,— 1,... — 1), and —p with the eigenvectors Vi, Wi 

(2 < i <n), given by 

Vi = Ci - Cj, 2 <i <n, 

Wi = e„+i - en+i, 2 <i <n. 

We thus see that (I22J) holds if p + g > and q> 0, holds if p + g < and 
q <0 and holds if q{p + q) <0. We thus obtain: 
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Theorem 10. Consider the system: 

2n 



61 = s(e,,p f{e,,0j) + q f{o^,e/)y i<i<n 

j = 1 j=n+l 
3 i 

n 2n 

ei = s[0,,qYf{^^^^J)+P E n+l<i<2n 



j=i j = „ + 1 

j 7^ i 

and assume |3|) holds. Define x by H2<^) . 

(i) If q{p + q) < and X 7^ ihe synchronized oscillation is unstable. 

(a) If p + g < and q < then the synchronized oscillation is stable if x < 

and stable ifx>^- 

(Hi) Ifp + q > and q > then the synchronized oscillation is stable if x > 
and unstable if x < 0- 
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